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Abstract. We consider a S-N-S Josephson junction modeled through the 
Ginzburg-Landau theory. When the normal material is sufBciently thin and 
the applied magnetic field is below the critical field of vortex nucleation, we 
prove to leading order that jump boundary conditions of the type predicted 
by de Gennes are satisfied across the junction. 

Introduction 

The superconducting proximity effect in a normal metal adjacent to a super- 
conductor has received a lot of attention by the physics community, see [TO] for 
a review of this phenomenon. This is also the setting of the Josephson tunneling 
effect for superconducting-normal-superconducting junctions (SNS), where a super- 
current flows through the normal layer provided that it is sufliciently thin. 
The physics literature contains several approaches to model the Josephson effect 
in the frame work of the Ginzburg-Landau theory of superconductivity. The first 
modeling in this context is perhaps due to the physicist de Gennes [12]. In the 
setting of [l^, the complex- valued wave function (whose modulus measures the 
density of superconducting electrons) and its derivative are related linearly on both 
sides of the normal material, in such a manner that the supercurrent is conserved 
through the junction. 

In this paper, we use a generalized Ginzburg-Landau energy functional presented in 
[9], which has proved to account rigorously to various physical aspects (c.f. [T9l[20] ). 
By working in the London singular limit (high k- regime), we justify asymptotically 
the modeling of [12] provided that the applied magnetic field is below the critical 
field of vortex nucleation, see Theorems 1 1 . II fc [ l . 51 

We also mention in this direction that another justification of the de Gennes mod- 
ehng is present in a paper of Rubinstein-Schatzman-Sternberg [23], who deal with 
geometric junctions (weak Hnks) in the framework of the Ginzburg-Landau theory. 
We hope to carry out in a forthcoming work a deeper analysis valid for higher ap- 
plied magnetic fields and which provides more details concerning the supercurent 
fiow and the distribution of vortices in the junction. 

1. Main results 

We move now to the mathematical set-up of the problem. Let fl — D{0, 1) 
denotes the unit disc in K^. Given R s]0, 1[ and i €]0, R[, we introduce the following 
partition of Q, 

n = SUN, 
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where 

(1.1) N ^{xen : dist{x,dD{0,R)) <£}, 

(1.2) Si = D{0,R~e), S2 ^ D{0,l)\D{0,R + i), 5 = 5i U S'a. 

We shall suppose that S is the cross section of a cylindrical superconductor with 
infinite height and that N is that of a normal material. By this way, we get a S-N-S 
Josephson junction. 

In Ginzburg-Landau theory [13], the superconducting properties are described by a 
complex valued wave function ip, called the 'order parameter', whose modulus \ip\'^ 
measures the density of the superconducting electron Cooper pairs (hence ^p = 
corresponds to a normal state), and a real vector field A = (^1,^2), called the 
'magnetic potential', such that the induced magnetic field in the sample corresponds 
to curl A. Since, the superconducting Cooper electron pairs can diffuse from the 
superconducting to the normal material in a normal-superconducting junction, we 
then have to consider pairs {ip, A) defined on il. 

The basic postulate in the Ginzburg-Landau theory is that the pair {ip, A) minimizes 
the Gibbs free energy, which has in our case the following dimensionless form [9] : 

(1.3) GcHi^J^A) - f ~ iA)^/j\' dx + ^ f {1 ~ dx 

+ / |Vpdx+ / \cMv\A- H\'^dx. 
^ Jn Jn 

Here, ^ = k is a characteristic of the superconducting material (filling S), H > 
is the intensity of the applied magnetic field and a > is related to the critical 
temperature of the material in TV. The positive sign of a means that we are above 
the critical temperature of the material filling TV. 
Minimization of the functional l|1.3p will take place in the space 

We will be interested in the analysis of the asymptotic behaviour of the minimizers 
of (|1.3p as e ^ (London Limit) and when the thickness of the ring TV is small by 
taking £ = £{e) < 1 as e 0. 

According to [15], the functional (|1.3p admits a minimizer {ip,A) in the space 
H. Our main result is the leading order asymptotic expansion of the jump of 
(V — iA)^p across the junction, i.e. across the boundary of TV, provided that the 
order parameter ip is not possessing vortices. 

In order to fix ideas, given a function / € H^{fl; C), we introduce the jump of f 
across TV by 

(1.4) [fU{e) = f{iR + i)e^')-f{iR-e)e^'), V0e[O,27r[. 

Our first result concerns the case of very thin rings, of thickness comparable with 

£. 

Theorem 1.1. Let Q = D{0,1), S and TV as in [Tl\) . [T^) . and (V'e,^e) be a 
minimizer of fl.3]) . Given d > and a > 0, there exists A > such that, if the 
applied magnetic field satisfies 

(1.5) F<A|ln£|, 
and if i = de, then IV'eI > and we have 

(1.6) lim 

(1.7) lim||[^,]A.||^,(g,^=0. 



'n{x) ■ (V - iAe)^^ 



_ exp(2yad) - 1 
— zya — r^— 

^ exp{2y/ad) + 1 



= 0, 
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Here, n{x) = - — r for all x eM.^ \ {0}, is the unit outward normal vector of any disc 

\x\ 

in R^. 



Remark 1.2. The regime concerning the applied magnetic field H in Theorem \l.l\ 
corresponds to that below the first critical field : When I11.10\} is satisfied, the order 
parameter ipe has no vortices in fl. 

On the other hand, it is well known (c.f. [2^^ that there exists A' > such that 
if H > A'jlnel, the order parameter ipe has vortices. However, we are not able to 
calculate the critical value Ac for which 

^^Ci Ac| lne| as e — > . 

This is due to the technical difficulty arising from the very rapid oscillations of the 
maximal superconducting density in N , see Section 2 for more details concerning 
this point. 



Remark 1.3. Notice that if one formally makes d ^ in lfl.6\) . one would obtain 
that the jump across N tends to 0. This agrees with experimental and theoretical 
predictions that the Josephson effect would be absent in junctions made up of very 
thin normal materials, see [9\. 

Remark 1.4. Let us introduce the vectors 

^- ^ , enix) ■ (V - lA,)^, { {R - ^)e^«) 



e 



nix) ■ (V-iA,)^, {{R + £)e'^) 



and the matrix 



1 2V^.^"P(2^^) 



Ma,d = I ^ cxp(2V^d) + 1 

1 

Then, ^.6\) and can be written in the equivalent form 

(1.8) lim \\X+ - MaMX; - 0, 

which justifies the boundary condition postulated by de Gennes in [12] . 
For thicker rings, we have a result analogous to that of Theorem ll.il 

Theorem 1.5. Let Q = £'(0,1), S and N as in [Tl\) . (EH), and {ipcA^) be a 

minimizer of fl.3\} . Assume in addition that I — £{e) satisfies, 

(1.9) e<^(e) (e->0), 3 c> 0, V e e]0, 1], ^(e) < celn | lne| . 

Then, given a > 0, there exists A > such that, if the applied magnetic field satisfies 

(1.10) H<Xe^p(-'^y^^)\\ns\, 



we have 



(1.11) lim 

£— ►O 



n{x) ■ (V - iA^)^^ 



2^^ 



= 0, 

L2(Si) 



N 

(1.12) lim||[^,]^||^,^g,j=0. 

X 

Here, n(x) — - — - for all x eM.^ \ {0}, is the unit outward normal vector of any disc 



in M^. 
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Notice that the result of Theorem 11.51 agrees with that of Theorem 11.11 when one 
takes formally d +oo in l|1.6p . 

We mention also that the result of Theorem ll.Sl is valid for magnetic fields slightly 
much lower than that of Theorem 11.11 Technically, this is due to the fact that we 
can not exclude the energy of vortices in TV (see Section 2). But heuristically, the 
reason is that the maximal superconducting density (i.e. the positive minimizer of 
(|1.3p for H = 0) in the regime of Theorem 11.51 is small inside the junction, hence 
the price of a vortex becomes for magnetic fields less than that of Theorem 11.11 
However, we were not able in this case to prove that the critical field for vortex 
nucleation has the order of exp ^—h/^lisl^ | lne|, though theoretical predictions in 
the physics literature say that vortices in N would be present for magnetic fields 
much below than that of a bulk superconductor, see p^. 




The result of Theorem 11.51 is still vaHd up to lengths £{e) — c*£|lne|, where 
c, G]0, 1[ is sufficiently small (this can be checked through minor modifications 
of the argument). However, since the corresponding magnetic field will be small 

<C 1, we do not focus on this last regim^. 

As a by-product of the analysis that we shall carry out, we get a result concerning 
the conservation of the current across N. The current is defined as the vector field 



(1.13) J, = (#s, (V - iA,)ij,) = {{i^,, id2 - lAl)^,,), (#„ (92 - iAl)^,)) , 



where (•, •) denotes the scalar product in C when identified with M^. 
The exact result concerning the current is the following. 

Theorem 1.6. In both regimes of Theorems ] 1 . 1\ & \T7^ the circulation of the current 
is almost conserved across the junction: 



where t{x) = ^ for all x eM.^\{0}. 

Finally, we comment on some past work concerning non-homogeneous supercon- 
ductors and pinning models. Unlike to our situation, pinning models previously 
considered correspond to a term in the Ginzburg-Landau functional having the 
form (a — with a being a smooth function. The first analytic work probably 

goes back to [I], where a — a^, assumed always positive, may depend on e with the 
restriction that it can not oscillate quicker than | lne| (|Vae| < C| Inej). Compared 
with our situation, the discontinuity of the coefficients leads to an order parameter 
oscillating between and Ug with iVuel > ^ in a thin neighborhood of dS. Later, 
in [6], the authors deal with the case when a has isolated zeros, and prove that 
vortices appear first at the zeros of a for magnetic fields having order 1. In [4], 
the authors allow the function a to have negative values, but the hypothesis of its 
smoothness forces the order parameter to be small on the boundary of the nor- 
mal side (hence, the surface superconducting sheath in the normal side is absent). 
Moreover, the situation in [1] is more related to the case of domains with holes and 
relies on the analysis carried out by the same authors in [3]. More recently, the 
author of the present paper showed in [18] that pinning of vortices is observed for 
magnetic fields near the first critical field when the function a is a step function. 
Let us also mention the very recent work of Alama, Bronsard and Sandier in [5] 



It is more likely that the regime H I is treated, without additional restrictions on the 
decay of H, by bifurcation arguments. In the present paper, we treat in detail the case H = 0. 




(1.14) 
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where a layered superconducting model has been investigated. There, the expres- 
sion of the critical field above which vortices are detected depends strongly on the 
thickness of the normal regions separating the superconducting layers. 

Organization of the paper. 

Section 2 is devoted to a description of the main points of the argument. 
Section [3] is devoted to a preliminary analysis of the minimizers of l|1.3p . 
Section |4] is devoted to an analysis of an auxiliary variational problem (this is the 
variational problem (|2.5p describing the Meissner state). 
Section [5] is devoted to prove a lower bound of the functional (|1.3p . 
In Section [6l we exhibit a vortex-less regime and we achieve the proofs of Theo- 
rems O and [Tl 

Finally, in Appendix [Aj we prove a uniqueness result concerning the solution of the 
canonical equation without magnetic field (this is Eq. (|2.ip ). and in Appendix [B| 
we discuss the difficulty behind the estimate of the energy of a configuration with 
vortices on the circle Sr. 



2. Sketch of proof 

Canonical Equation in for the case without magnetic field. 

A general technique to tackle asymptotic problems, (successfully applied by Helffer 
& co-authors for linear problems (cf. e.g. [I6j), and Lu-Pan for non-linear problems 
(cf. e.g. p2])), is to understand the canonical situation. So we consider the case 
without magnetic field, H — 0, and work with the particular domains, 

N = Rx]- £Jl S = R^\N. 

When £{e) ~ de, this leads to the equations: 

r-Au=(l-u2)u in R X {R\ 

(2.1) < -Au + au= inMx]-d,d[, 

[ dx2u{-,t±) = dx^u{-,t^), u(-,i±) = u(-,tzp), t = ±d, 

Since Eq. (|2.ip arises as a limiting equation of a critical point of (G-L) , we focus 
on solutions in the class 

C ^ {uG Hl^{R^)nL°°{R'^) : u>0}. 

Using the argument of [19^ Section 4], Eq. I|2.ip admits in C a unique solution 

]R2 9 {xi,X2) ^ U{X2) 

given by 

r f3{a,d) exp(72|a;2|)-l , , > ^ 

(2.2) Uix2)^l /3(a,d)exp(V2|x2|) + l ' ' 

[ ^(a, d) [exp(^|a;2|) +cxp(-^|a;2|)] if |a;2| < d. 

where the constants (3{a,d) > exp(— \/2d) and A{a,d) > are determined explic- 
itly, but the only important remark is that 

U'(d) U'(-d) „ _exp(2V^d) - 1 
— — 2va- 



U{d) U{~d) ^ exp(2V^d) 4- 1 ' 
hence the conditions of de Gennes are verified 

(2.3) )..r[o exp(2V^d) + l I 
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The case without magnetic field in a bounded domain. 

Now we return to minimizing (|1.3p when H = 0. We prove in |191 Theorem 1.1] 
that (|1.3p has, up to a gauge transformation, a unique minimizer (we,0) where 
Me £ H^{n) is a real- valued function, < Ue < 1 in (for e small enough), and 
solves the equation 

1 2\ 

Then, by a blow-up argument, we generaHze l|2.3p asymptotically as e ^ 0, see 
Theorem [ 



The case with magnetic field: A vortex-less regime. 

Now we return to minimizers {-ip, A) of (|1.3p when H > 0. Following an idea of 
Lassoued-Mironescu [21], we introduce a normaHzed density 



Then, \ip\ < 1 and we are led to work with the functional (see Theorem I3.5P : 

+ / l'^n'da;+ / \cnTlA-H\^dx. 

^£ Js Jn 

Then, following |25[ [27] . we construct a family of disjoint balls {B{ai,ri))i (see 

Proposition I5.2p such that X^i ''i — |lii£|~^° and 

{xen : |v?(x)| < 1- |ln£|-2} c |jB(a„r,). 

i 

This permits to obtain, for a given number a g]0, a lower bound of the energy 
(see Theorem I5.4p : 

(2.4) Te,Hi^,A) > Moie, H) +^{2nami{e)\\iie\ - 2H)di - Cff|lne|- 

i 

where C > is an expHcit constant, di is the degree oi •^/\ip\ on dB{ai,ri), 



-10 



mi{e) = min ^^(x) , 

x£B(ai,ri) 



and 



(2.5) Afo(s,i7)^^Jnf^^^^^^..(M). 

By comparing i|2.4p with what is existing in the literature (c.f. |26j), we suspect 
that the lower bound l|2.4p is not optimal in the sense that a should be equal to 
1. This restriction is actually due to the particular expression of T^^h, where a 
penaHzation term for \ is absent in the energy of N . 

The infimum in l|2.5p is achieved by a unique vector field A ~ ^V-^h^, where 
: fl I— >]0, 1[ satisfies a London equation with weight (see (|4.4p ). Thus, we get 
the upper bound 

J'e,H{'P,A) < Moie,H). 

When this upper bound is matched with l|2.4p . we deduce that all the d[s are equal 
to provided that 

H < an ^inf rni(e)^ | lne| . 
If £{e) ^ e and \ai\ = R, we have by Theorem 13.41 that 

mi(e) ~ exp ' 
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Since we can not exclude the possibility of a vortex ball centered on the circle §)j, 
we restrict ourselves when £{e) £ to magnetic fields H satisfying 

if<Aexp(-^^)|lne| 

in order to insure the absence of vortices. 

Now, in the absence of vortices we get an energy estimate (see Theorem 16. ip 
/ (|(V-iA)(p|2 + |curM-i7/i,|2)dx+4 /(1-I'^n'dx<l (e 0) . 

Then, we implement L^-estimates for the equations of ip in order to deduce that 

Mtv^O inL2(S^), n{x) ■ (y - iA)ip ^ mL^{dN), 

which permits us to deduce Theorems 11.11 and 11.51 see Section 16.11 



3. Preliminary analysis of minimizers 

3.1. The case without applied magnetic field. This section is devoted to a 
summary of the main results obtained in [19] which deal with minimizers of (|1.3p 
when the applied magnetic field H = 0. 

We keep the notation introduced in Section 1. Upon taking A ^ and H ^ in 
(|1.3p . one is led to introduce the functional 



(3.1) Geiu):^ [ \\/u\^dx + ^ [ {l-u^)^dx + ^ [ u^dx, 

Jn ^£ Js ^ Jn 

defined for functions in If^(r2;M). 
We introduce 

(3.2) Co(e)= inf g,{u) . 

The next theorem is a summary of Theorem 1.1 in |19) . 

Theorem 3.1. Given a > and d > 0, there exists Eq such that for all e e]0,eo[, 
the functional 113.1]} admits in H^{n;M.) a minimizer G C^(S') UC^(iV) such that 

< Ue < 1 in n. 

Furthermore, with our choice of the domains Q,, N and S in (GUP and (HE), the 
function Ue is radial. 

If H = 0, minimizers of fl.3]} are gauge equivalent to the state {ue,0). 

Let us just mention why we focus only on the regime e ^ in the statement 
of Theorem 13.11 Notice that u = is a critical point of the functional l|3.ip , so we 
would like to exclude the possibility that this critical point is stable. To this end, 
we define the following eigenvalue : 

Ai(a,d,£) = inf 1^ (^|V(/.|2 - dx 

(|V</.|2 + 4|</.|2) da; : G H^n), Uh^^^^ = l] . 

Then when Ai(a, d, e) < 0, the corresponding eigenfunction of Ai(a, d, e) provides 
us with a test configuration whose energy is below that of m = 0. On the other 
hand, this last condition of the sign of Ai(a,TO,e) is easily verified when e ^ 0, 
thanks in particular to the min-max principle. 

Let us recall the notation of the jump across N introduced in (|1.4p . Using a blow- 
up argument and a result concerning uniqueness of critical points of the functional 



7. 
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(|3.ip in the model case of the entire plane, we establish now Theorem 11.11 in the 
case when H — 0. 



Theorem 3.2. Let he the positive minimizer of IfS.l]) introduced in Theorem \3.1 

Then, if i{e) - de, we have 



(3.3) 
(3.4) 



lim 



n{x) ■ VWe 



■ exj>{2^/a d) 



^ exp(2Vaa) + l 



0, 



L~(Si) 



lim II [us]jv|L^(si) =0, 



(3.5) 



lim I sup Ue{{R± £)e*^) - A(a, d) 



0. 



Rere, n(x) = - — - for all a; G \ {0}, and A(a, d) > is an explicit constant. 

\x\ 

On the other hand, if l{s) ^ e (this covers the regime (1.9]} }. then we have, 

n(x) ■ Vu^ 



(3.6) 



lim 

e-+0 



N 



0, 



L~(Si) 



(3.7) 



lim sup \uJ{R± t)e''^) ~ A{a) | = , 



where A{a) > is an explicit constant. 

Proof. Let us treat the case when £{e) = de, d > 0. Let (r, 6) be polar coordinates, 
< r < 1, —TT < 9 < TT, and set 

t = R-r, s = R9. 

Given sq G [— -Rtt, i?7r[, we define the rescaled function, 

u,{s, t) = u, ((i? - et)e''^'-'">^/^^ , 
The equation of Ue becomes: 



R-1 1-R R R 

< t < , — TT— < S — So < TT—. 



-AeUe + amue = 0, ^^<t<0, |s-so|<7rj, 
^{■,t±) = ^(''^t)' Uei-,t±) = u(-,%) for t = ±d. 



where 



A, = 1 - e 



R 



(R-et) 

Now, by elliptic estimates, the function converges to a function u in W^^^ {R'^) . 
Furthermore, u solves (|2.ip in C, and by [T9i Lemma 5.2], there exist constants 
ko, Co > such that u{0, ko) > cq. Thus, we conclude by Theorem lA.ll that u{s, t) = 
U{t), where U is given in l|2.3p . and therefore, by coming back to the initial scale, 

Ue(s, t) — U 



Vfce {0,1,2}, limf 



0, 



lV'=.==({|s-so|<7rfle, |t|<(l--R)e}) 



and the convergence is uniform with repect to sq G [^tt-R, vri?[. This yields p.3p 
(13:511. 



The expression of A{a, d) is given explicitly in the Appendix. 
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The statements concerning the case when £{e) ^ e follows from [T9| (5.20)]. □ 



We shall need the following remarkable properties of u^, that distinguish the 
different regimes considered in this paper. 

Lemma 3.3. With the notations and hypotheses of Theorem \3.S\ if £ ~ de, there 
exists an explicit constant c(a, d) > such that 

Ue{x) > c{a, d), \f X Q Q, . 

For the case when £{e) satisfies (|1.9p . the asymptotic behavior of becomes 
completely different in the sense that it is close to inside N. 
In order to be precise we introduce the function : 

/3exp(V2i) - 1 



(3.8) 



Vit)^ 



{t>0), V{t)=Acxp{y/Et) (i<0), 



/3exp(\/2i) + 1 
together with the 'signed distance' to the boundary of S, 

(3.9) tsix) =distix,dS) (x € S), tsix) = -distix,dS) {x e D{0,1)\ S). 
Here the constants (3 and A are given by : 



(3.10) 



Theorem 3.4. Assume that £{e) satisfies il.9\) . Then, we have 

(ts{x) 



(3.11) 



V e 



= o(l) (e^O), 

where the functions ts and V have been introduced in f3.8\) - (3^} . 
Moreover, there exist a positive constant Eq and a function ]0, 1] 9 e i-^ g{e) e]0, 1[ 
such that g{e) <C 1 and for all e g]0,£o], one has the estimate 
(3.12) 



{A-g{e))e-xp 



\/ats{x) 



< u,{x) < {A + g{s))exp 



^ y/ats{x) 



y X e N. 



Proof. The asymptotic behavior in l|3.1ip has been obtained in [19]. We have only 
to prove the improved estimate in N, i.e. I|3.12p . 

Let us show how one can obtain the lower bound. Let us introduce the function : 

Ve [x) — C exp I 

where C and S are positive constants to be specified later. 
Let us recall that by definition, the function ts is written as 

" i? - ^(e) - l^l if R ~ £{e) < \x\ < R, 
\x\^R-£{e) if R<\x\<R + i{e), 

where the constant R e]0, 1[ has been introduced for defining S and N, see (|1.1 
and ifO]) . 

Consequently, the function ts is smooth in each of the following two annuh of N: 

N^={xe N : R- £{e) < |a;| < i?}, N+ ^ {x e N : R < \x\ < R + £(e)}. 
One then checks easily that 



. , .a S 

-A{Ue-Ve) + ^{Ue-Ve) = ^ 



1 " _L ^1 I 



1 



52 d{R-£is)) 



in N- 



±- 
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It is a result of the asymptotic formula p. lip that there exist a constant Eq and a 
function ]0,£o] 9 £ i-^ /(e) G]0, 1[ such that /(e) <C 1 as e ^ and 

\ue\g„ - A\ < f{e), ee]0,eo]. 

Therefore, gathering all the above remarks, we get for 



and when e £]0, eo]. 



a 



(3^3^ I -A{u,-v,) + —{u,-v,)>0 mN± 

[ Ue{x) — Ve{x) > On ON . 

Two cases may occur regarding the gradient of on the circle \x\ = R, either 
<(i?) < or u'^{R) > 0. 

If <(i?) < 0, then we get in addition to ([313]) 

-^{ue-v,)>0 oii{dN+)nN. 

Here, we recall that vn± denote the unit outward normal vectors of the boundaries 
of 7V±. 

Therefore, the application of the strong maximum principle yields that 

— > in iV+ . 
This last lower bound when combined with l|3.13p yields 

I Uj(x) — v^{x) > on 3iV_ . 

Hence by the strong maximum principle, Ue — > in iV_. Therefore, we deduce 
that 

— We > in , 

which is nothing but the lower bound we wish to prove for the function u^. The 

same argument holds when u'^{R) > 0, but by changing the roles of and 

The proof of the upper bound follows the same lines above, so we omit the details. □ 



3.2. The case with magnetic field. This section is devoted to a preliminary 
analysis of the minimizers of i|1.3p when H ^ 0. The main point that we shall show 
is how to extract the singular term Co(e) (cf. I|3.2p ) from the energy of a minimizer. 

Notice that the existence of minimizers is standard starting from a minimizing 
sequence (cf. e.g. |l5j). A standard choice of gauge permits one to assume that 
the magnetic potential satisfies 

(3.14) divA = inn, v-A = Q on dn, 

where v is the outward unit normal vector of dVl. 

With this choice of gauge, one is able to prove (when the boundaries of Vl and TV 
are smooth) that a minimizer (V', ^) is in C^(ri;C) x C^(ri;R^). One has also the 
following regularity (cf. |19[ Appendix A]), 

V-e C2(^;C)U (7^(77; C), A C'^(S\R^)lJ C^(N]R^). 

The next lemma is inspired from the work of Lassoued-Mironescu (cf. |21j). 
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Lemma 3.5. Let (ip.,A) he a minimizer of il.3\} . Then < IV'I < in fl, where 
Ue is the positive minimizer of ^3.1]) . 



Moreover, putting if = then the energy functional fl.3\} splits in the form : 

(3.15) g,^H{^,A) ^ Coie) + J^e,Hi'P,A), 

where Co(e) has been introduced in ^3.2]) and the new functional T^m is defined 
by : 

(3.16) T,,h{v.A)^ [ ul\iV ~ iA)^\^ dx 

+ ^ ( u^(l- |^|2)2dx+ / \cni\A~H\''dx. 

Proof. 

The equality p.l6p results from a direct but some how long calculation, which 
permits to deduce in particular that is a solution of the equation 

-(V ~ iA)ul{V ~ lA)^ = ls^(l ~ Iv^nV • 

Proof of IV'I < u^. 

It is sufHcient to prove that \ip\ < 1. We shall invoke an energy argument which we 

take from |11| . 

Let us introduce the set 

n+^{xen : \ip{x)\ > 1}, 
together with the functions (defined in : 

/ = ^=[1^1-1]+/- 
m 

Then, it results from a direct calculation together with the weak-formulation of the 
equation satisfied by Lp that 



= / (|V|^|p + (|^|-l)|^||(V-z^)/p)u?dx 

+1 / {i+\^\){i-\^\r\^\)utdx. 

^ Jn+ns 

Therefore, this yields that C N and that |V|(p| | = in Hence, |(^| is con- 
stant in each connected component of which shows that \ip\ = 1 in This 
contradicts the definition of n+ unless n+ =^. □ 

The estimate of the next lemma is very useful for exhibiting a vortex-less regime 
for minimizers of (11.31). 



Lemma 3.6. Let {4>, A) be a minimizer of (pljj. There exist constants C > and 
Eq s]0, 1] such that, if the applied magnetic field satisfies <C -j, then we have 

\{W-tA)ij\<j, yee]Q,eo]. 

Proof. The proof is essentially due to Bethuel- Riviere [8], but we include the main 
steps for the reader's convenience. 

Since |(V — i^)?/;! is a gauge invariant quantity, we assume that we are in the 
Coulomb gauge (|3.14p . 

Let us assume that the conclusion of the lemma were false. Then there exists a 
subsequence, denoted again {e} and points (Xg) £ fi such that 

(3.17) e\{V - iA)'il!{xe)\ ^ oo . 
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We define the rescaled functions 

Ve{x) = ip{x^ + £x), B^{x) — £A{x^ + ex) , 
together with the rescaled domain 

= {n-x,)/e. 

Notice that {v^,B^) satisfies the equations 

-Ave + 2iB,-W = [{l-\v,\^)ls-alN + \Be\'']v, in 17^, 

— ASg = {iv^ , Vuj — iB^Vg) in fig, 

curlSj = e^H^ on dfl^, 

n ■ Vfe = on dH.^ . 

Notice that is a weak solution in fl of the first equation above because of the 
transmission conditions 

where 

are the 'interior' and 'exterior' trace operators. 

With the choice of gauge in p.l4p . we get by Poincare's Lemma that ||^||ifi(o) < 
C|jcurl A||2,2(Q'). Then, by elliptic estimates and the assumption H ^ ^, we 
have 

\\A,\\H2(^n) < C||curlA||L2(o) < -. 

Hence, by the Sobolev embedding theorem, we have in the new scale, 

lim ll^elU-(n,) = 0. 

Now, by Lemma [331 \ve\ < < 1, hence, the right hand side of the equation of 
becomes bounded. Therefore, 

\\Av,\\L,^Bn) <Cr + 2||B,||ioo(o.) X ||Vfe||LJ.(Bfl), V p > 2 , 

where B^ is any fixed ball of radius R. 

By elliptic regularity theory, becomes bounded in W^'P^Bh) for all p > 2, hence, 
by the Sobolev embedding theorem, in C^'°'{Bpi) for any a €]0, 1[. Since C^'"(i?ij) 
is compactly embedded in C^{Bji), we get by a diagonal argument the existence of 
a function v e C^(M^) such that, upon extraction of a subsequence, converges 
to V locally in . 

Now, we know from the equation of B^ that ABg is locally bounded in L°°. So 
again, by elliptic estimates, and since ||i?e||Loo as e — + 0, we get upon extraction 
of a subsequence that B^ converges to locally in . 
Therefore, we get, by returning to the initial scale, 

e\{V - iA)ij{x,)\ = \{V - iB,)v,m 
is convergent, hence contradicting l|3.17p . □ 

Now, Lemma \3M permits to conclude the following result. 

Lemma 3.7. Assume that {ip,A) is a minimizer of ( f and let tp = There 
exists a constant /xq > such that if 

4 / {l-H^fdx<^io, 
£ Js 

then \(p\ > in S. 
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Proof. Lemma \3M and the diamagnetic inequality together yield that 



Now, since 



we deduce that 



IVIV'I I < |(V-i^)V'| < -, infl. 

e 



\Vu,\ < - 
e 



\V\ip\\<- \nS. 

£ 



Thus, the result of the lemma becomes a consequence of [TJ Theorem III. 3]. □ 



4. Analysis of the Meissner state 

Let us recall the definition oiu^ and Co(e) in Theorem 13. II and l|3.2p respectively. 
This section is devoted to the analysis of the following variational problem : 

(4.1) Mo(s,i/)^^^min^^^^,.(..,^^ 

Since the function is real-valued, one gets, for any vector field A, the following 
decomposition : 

g,,H{u,,A)^Co{e)+ f {\Au,\^ + \cmlA~H\^) dx. 
Jn 

Putting further 



A = HA, 



(4.2) Jo(£)= inf 

we get that 



/ (|^MeP + |curM-lp) dx 
Jn 



Mo{e, H) = inf ^e.ff (we, A) = Co{e) + H^Me), 



and we are reduced to the analysis of the variational problem l|4.2p . 

Starting from a minimizing sequence (cf. [26]), it is standard to prove that a 

minimizer A^ of l|4.2p exists and satisfies the Coulomb gauge condition: 

div = in ri, n ■ A^ = {] on dVl, 

where n is the unit outward normal vector of the boundary of fi. 
Notice also that A^ satisfies the Euler-Lagrange equations : 

(4.3) V-^ curiae = in 17, curM^ ==1 on dn. 
Here V"*^ = {—dx^Tdxi) is the Hodge gradient. 

Putting — curiae, we get from the first equation in (|4.3p that Ag, — ^^V^/i^. 
We get also that satisfies the equation: 

(4.4) - div (^\^h^ + K = in Vl, h, = I on dn. 

Lemma 4.1. The function satisfies < < I in fl, and it is the only function 
solving |73P- 

Moreover, there exist constants co,eo g]0, 1[ such that, 

(4.5) Co < - l||L-(f2) < 1, V£e]0,£o]- 
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Proof. The property that < hg < 1 and the uniqueness of /i^ are direct appHca- 
tions of the Strong Maximum Principle. 

Let us now prove l|4.5p . Assume by contradiction that there exists a sequence 
converging to 0, still denoted by e, such that 

(4.6) lim||/i,-l|U^(n) =0. 

Let us take a compact subset K C S (independent of e). Due to the asymptotic 
behaviour of (it remains exponentially close to 1 in K, see Theorem [SH] and [lH 
Proposition 5.1]), it results from l|4.4p that is bounded in the C^-norm of K. 
Thus, one can extract a subsequence of h^, still denoted by he, that converges to a 
function h G C'^{K). The function h satisfies the limiting equation, 

-A/i + /i = in K. 

Coming back to (|4.6p . we get that h=l'mK, hence not a solution of the limiting 
equation. Therefore, the assertion in l|4.5p holds. □ 

The next results concern the case of our particular domains in (|l.ip and (|1.2p . 

Lemma 4.2. With the assumptions II 1.1]) and ^.2]) . the function h^ is radial, i.e. 
he{x) — he{\x\), with he being an increasing function. 

Proof. That h^ is radial follows by the uniqueness of the solution of (|4.4p and by 
the fact that is also radial. 
The solution he being radial, i.e. 

he(x) — he(\x\), V a; e il, 

let us show that the function he is increasing. For simpHcity of notation, we shall 
remove the tilde and write he for he. Notice that he satisfies the differential equa- 
tion : 

(4 7) I -K{r) - \h',{r) + 2^ K{r) + ul{r) heir) =0, < r < 1, 

h' (r) 

Let us calculate h'^{0). Since h'^{0) — 0, we have /i"(0) = 1™ ~^ — • Substituting 
in l(4?7l) . we get that 

(4.8) K{0)^^u'emheiO)>0. 

Let us introduce the even extension of he, namely the function 

heir) (r > 0), 



Mr) = 



hei-r) (r < 0). 



Then fe satisfies the equation, 

(4.9) - f:ir) ~ -f'Ar) + 2^ /:(r) + u?(r) A(r) =0, re] - r^, r2[\{0}, 

I'l Ue\r) 

and it attains a local minimum at 0. We emphasize also here that Ue denotes the 
even extension of Ue . 

If ro e] — 1, 1[ (with ro 7^ 0) is a critical point of fe, then it follows from (|4.9p that : 

/;(ro) = u2(ro)/e(ro)>0. 

If ro = 0, the conclusion /"(O) > still holds, thanks to ijTSl) . 
Now these observations lead to the conclusion that fe attains its minimum at a 
unique point, and that this point is the only critical point for fe. As we know that 
/^(O) = 0, we get that fe attains its minimum at and that it is increasing in [0, 1[. 
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This achieves the proof of the lemma. 



□ 



The next lemma plays a distinguished role in the control of the minimizing energy 
of 'vortex balls'. 

Lemma 4.3. There exist constants C > and eo > such that 



(4.10) 



ui 



< 1, VeG]0,eo]- 

L°°(n) 

Proof. Notice that by Lemma \A.2\ is radial. Then the equation for he can be 
written in the form: 

1 h 



^ M + /i,(r) = 0, Vre]0,l[ 



Integrating this equation between and r e]0, 1[ and using the fact that is 
increasing, h'^ > 0, we obtain: 

''l) (r) < J\e{7)d¥< r||/i,|U^([o,i]) < 1, 
which is the result of the lemma. □ 



5. Lower bound of the energy 

5.1. Construction of vortex-balls. We borrow some notation used in [26]. For 
a set [/ C we denote by r{U) the radius of U, that is the infimum over all finite 
coverings of U by open balls Bi, B2, ■ ■ ■ , Bk of the sum ri + r2 + • ■ • + . The 
important property is that: 

r{U) < \H\dU), 
where is the one-dimensional Hausdorff measure. 



From now on, we shall always work under the following hypothesis: 

(H) 3c> 0, V e e]0,l], £(e) < cln|lne|, H<cme\lne\ 

where 



(5.1) 



TOe = in£ug(x) 



Notice that when £{e) = de, and ^ are bounded. When i{e) ^ e, we have 
by Theorem [331 

^^^^(^)\l + o(l)), (S-.0), 



rrie ~ A{a) exp ^— 
where A{a) > is an explicit constant. 

Lemma 5.1. Let (V', ^) 6e a minimizer of S1.3\) and Lp — — . Then, under the 
hypotheses (H) above, there exist constants C > and £0 > such that, for all 



{\\ne\ + |lnme|) 



^,1 

£0 



we have 
(5.2) 



r[{xen : \ip{x)\ <l-5} ]<C 



and e s]0, Eq], 
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Proof. We have the following decomposition of the energy, 

g,.H{^,A) = Co{e)+T,,H{v,A), 
where the functional J^e,// >1) has been introduced in p.l6p . 
Using {u^, ^V"'"/ie^ as a test configuration for the functional (|1.3p . we get 

(5.3) T,^H{'il^,A)<ZH''. 
We infer from il5.3L 



/ |(V-a)(^pds + -i^ / {\-\^\'' fdx<2cm-^H\ 
Jn ^£ Js 



where is introduced in (|5.1 
By Lemma [3T5I \(p\ < 1, hence 



(5.4) 



/ {l-\^\^fdx < 2\N\ < 4tt£{s) < Ce\\nm,\. 
Jn 



Now, by l|5.4p and the diamagnetic inequality, |(V — iA)ip\ > |V|(,5| |, we deduce 



(5.5) 
where 



^ p dx + ^ ^(1 - |^|2)2 dx<M, 



e = V£, M ^2c{mJ^H^ + \lnm,\) . 

Since ip is a i/^-function, then it can be approximated in the L°°-norm by means 
of smooth functions. This permits us to conclude l|5.2p from Proposition lC.il □ 



The next proposition provides us, as in |26j, with the construction of suitable 
'vortex-balls'. 



Then, under 



Proposition 5.2. Let {tp,A) be a minimizer of 111.3]} and ip 

the hypotheses (H), for each p e]l,2[, a g]0, l/2[, n G N, there exist constants 
C > 0, 7 s]0, l/2[, and for each ^ ?y ^ 1, there exists a family of disjoint balls 
{B{{ai,ri)}i satisfying the following properties: 

(1) w = {xen : |(^(a;)| < 1- |lne|-"} C U,B(a„rO. 

(2) J2r,<V- 



(3) Letting di be the degree of the function p/\p\ restricted to dB{ai,ri) if 
B{ai,ri) C fl and di — otherwise, then we have: 



(5.6) 



ul\{V - iA)ip\^dx + 



B{ai,ri)\u 



\c\xr\A- dx > 



27r|d,- 



B{ai,ri) 



mm u. 



(4) 



277 ^ diSai - curl(A + {iip, V ap)) 



In — - C|lnf 



< max(|lne|2-2»,,yT). 



Proof. Let us take (5 = |lne| " in the statement of Lemma [5.11 We emphasize 
that under the hypothesis (H), our choice of 5 verifies the hypothesis of Lemma [5?T| 

(|ln£| + |lnm,|)./5,l 



namely, 5 E 

Points (1) and (2) of Proposition 15.21 are now direct consequences of Lemma [5.11 
and the first point of Proposition IC.2I 

Let us prove now Point (3). By the estimate on r{w), we get for a given a €]0, 



JUMP CONDITIONS FOR JOSEPHSON JUNCTIONS 



17 



r{w) < provided that e is sufficiently small. 
Now, notice that in B{ai,ri) \ w, 

|(V-zAVp>|^|2 

Hence, we get the desired conclusion by applying the third point of Proposition [C2] 
to the function ^/\(p\. 

Notice that the function ip and the balls B(ai,ri) satisfy the hypotheses of Propo- 
sition [CH] with M = ©dlnep). The application of this proposition yields the 
conclusion in the last point of Proposition 15.21 □ 



(V - iA) 



> -2|lne| 



Remark 5.3. When performing the previous argument with 5 = and [3 sjO, |[ 
sufficiently small, we get improved remainders in (3)-(4) of Proposition \57^ a power 
of e, hut valid for larger values ofrj. This permits to treat the case when the thickness 
^(e) < ce| lne| and c e]0, 1[ is sufficiently small. 

We follow the usual terminology and call the balls constructed in Proposition [SH] 
'vortex-balls'. 



We conclude with the following theorem. 



Theorem 5.4. Let {ip,A) be a minimizer of (LSp and ip = — . Then, under the 

hypothesis (TV), for each a e]0, l/2[ and n e N, there exist a constant C > and a 
family of disjoint balls {B{{ai,ri)}i such that : 

(1) («'G]o,i-a[); 

i 

(2) \(p\ > \ onn\(J.,B{a„r,). 

(3) Letting dt be the degree of the function yy/\ip\ restricted to dB{ai,ri) if 
B{ai,ri) C fl and di = otherwise, then we have: 



(5.7) 



di>0 



mm u, 

3(a.,n) 



llnf 



2H 



d, - CHI In el 



(4) Y.^\d^\<Cm,\\ne\ 



Proof. Applying Proposition 15.21 with rj = e" and with a replaced by /3 e]0, ^[ to 
be chosen sufficiently close to ^ , we get a family of balls satisfying in particular the 
first two assertions of the theorem. 



Total degree. 

We start by proving an upper bound on the total degree \di\. The starting point 
is by noticing that there exists a constant c > such that 

T,,Hiip,A)<cH\ 

Then, by applying Point (3) of Proposition l5.2l with 77 = e" , we deduce the existence 
of a constant c > such that 

27rme|lne|^ < cH'^, 

i 

hence Point (4) of the theorem is proved. 



AYMAN KACHMAR 



A rough lower bound of the energy. 
We put 

B{ai,ri)Cfl 

H 

12 



JLJ- 



Since \lp\ < 1, then 



and consequently, we have 



Te.Hiip.A,n) > {ul\j\^ + IcurM- dx. 
Jn 

Now a direct calculation yields, 

T,.H{f,A,n) >H^Jo{e,n) + 2H [ {h,-l)[cm\{A'+j')] dx 

- 2H 



f [ih,-l)cUTlA' -f -V^he] dx, 

JuiB{ai,ri) 

1 



where 



l/i. - IM dx 



SI V" 



Using Lemma 14731 and the Cauchy-Schwarz inequality, it is easy to prove that 



[ [{he - 1) curM' - / • V^/ie] dx 



and 



Jo{e)-Jo{e,n) <CY,n- 

i 

Therefore, we obtain 

(5.8) J^e Hi'f, A, h) > H^Jo{e) + 2H [ {h^ - 1) [curl(A' + j')] dx - CH^ V r,. 

Jn 

Since curl(A' + /) — curl {A + {iLp, "S/ av)), then by Point (4) in Proposition [521 we 
rewrite the above lower bound in the form 

(5.9) T,^H{ip,A,n) > H^J^{e)+A^THY,d^{he - l){a,) - CH\\ne\-'^. 

i 

We have also by Point (3) of Proposition 15.21 Va e]0, W and when e is sufficiently 
small. 



(5.10) ^ T,^n{^,A,B{a,,r,))>2T^oi^\d,\\ min w^llnel. 
Therefore, we obtain the lower bound 

(5.11) TeuWA) > -ffVo(e) + 27ra|ln£|y ( min u^j 

i \Bia^.ri) ) 

+ATrH^di{he - l){a,) ~ Ci/llnep". 

i 

Since < h^ < 1, the above lower bound is sufficient to deduce (15.71). 



□ 
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6. Proofs of main results 

6.1. A vortex-less regime. Let us recall the definition of the constant rrie intro- 
duced in (|5.ip . We recall also that A) always denotes a minimizer of (|1.3p and 
that ^. 

The aim of this subsection is to prove the following theorem. 

Theorem 6.1. There exists a constant A > such that i/ £ < ©(e In | lne|) and if 
the magnetic field satisfies 

H < ATOe|lne| , 

then 

Ivl > ^ ills', 
and we have the energy estimate as e ^ 

f ^tA')ip\^ + \cm\A- Hh,\'') + 4 / (1 - It^H' < 
Jn £ Js 

Here 

A' ^A-^rV^h, 

and hi; is the function introduced in jji-j^ . 

One essential step towards the proof of Theorem 16.11 is a further useful splitting 
of the energy due to Bethuel- Riviere (cf. [8j). 

Lemma 6.2. Consider {u,A) e H^{n;C) x i/i(f7;R2) and define 

A' = A^^V^h„ 

where and /i^ are introduced in Theorem \ 3.1\ and Hj-j^ respectively. Then we 
have the decomposition of the energy, 

J'e.Hiu^A) =iJVo(e)+ / {ul\{V ~iA')u\^ + \cm:\A'\'^)dx 
Jn 

+ \ [ util - \u\Ydx + 2H [ (h, - 1) 
£ Js Jn 

+ [ ^{\u\^-l)\Vh,\^dx. 
Jn 

Here, the functional J-^^h and the energy Jo(e) are introduced in f3.16\} and {4-^ 
respectively. 



dx 
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Proof of Theorem \6.1[ 
Existence of X . 

Let us choose A > in such a way that when the magnetic field satisfies 

H < AmJlnel , 



all the degrees di given by Theorem 15.41 are equal to zero, so that the energy of a 
minimizer becomes close to that of the Meissner state. 
Matching the lower bound l|5.7p with the upper bound 

(6.1) T,Miv,A)<H^Me), 
we obtain, 

[am,\\ne\-2H]d^ + CHY^ - Ci?| lne|-" <0. 

di>0 di<0 

By our choice of H and i{e), we get that for any n g N, there exists a constant 
C > such that 

[am,\lne\ - 2H]di + CH Y < Uner". 

di>0 di<0 

Now we get by choosing A < 2a that, for all A < Ai, 

0<i/min(^-2,c)^|d,| < llnep". 

i 

Therefore, for all i, di — and there does not exist vortices. 
Moreover, the lower bound (|5.7p becomes when A < Ai, 

(6.2) T,,H{^,A)>H^Jo{e)~C\\ne\-^, 

hence when combined with the upper bound (|6.ip together with the energy expan- 
sion of Lemma lG^ and Item (4) of Proposition l5.2l we are able to deduce the energy 
estimate of Theorem 16.11 

Proof of \ip\ > i. 

Now, as we have the energy estimate, we deduce that 

1 

72 



(l-|(^ndx«l (e^O), 



5 

where we have also used that, in S, > ci for some explicit constant ci > 0. 
Therefore, we infer from Lemma [3771 that \Lp\>^m.S. □ 

It results from Theorem 16. II a uniform estimate of 

Corollary 6.3. With the notations and hypothesis of Theorem \6.1\ and if the 
magnetic field satisfies 

H < Am^l lne| , 

then we have as e — > 0, 

Proof. Notice that since and are radial, then 

Vu^{x) = n{x)u'^{\x\), Vhs{x) = n{x)h'^{\x\), 
where n{x) = |f| for all x € \ {q}. 

Therefore, the vector field A' satisfies the properties inferred from A, 

div A' = in n, n{x) ■ A' = on dn . 
Thus, it is a result of the Poincare Lemma that 

\\A'\\HHn) < C||curlA'|U.(o) 
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for some explicit geometric constant C > 0. 

Now, from the energy estimate of Theorem l6.ll we conclude the result of the corol- 
lary. □ 

The next lemma is now an essential step in proving Theorems ll.im.5l 



Lemma 6.4. With the notations and hypothesis of Theorem \6 . 1\ and if the mag- 
netic field satisfies 

H < Xm^ \ lne| , 

then we have as e ^ 0, 



and 



s\\n{x) ■ {V -iA')ip\\^^^gj^^ < 1, 
II [<^W|L2(si^) « 1- 

Here n{x) = - — - for all a: G \ {0}^ and [•]n is the jump across N introduced in 



Proof. Let us introduce the two domains 

Si = D{0, R-£)\ D{0, R/2 - £), S2 = I?(0, [R + l]/2 + ^) \ D{0, R + i). 

Notice that the domains Si and ^2 can be identified by those corresponding to 
1 = via the translations 

T± -.x^Re'"^ {R±e)e''^ . 

In order to prove the first statement of the lemma, we have only to establish (thanks 
to the trace theorem), 
2 

(6.3) 5^£|||n(x).(V-*A')^|||^.(5^)«l (e^O). 
1=1 

Notice that, by the energy estimate of Theorem 16.11 we have only to estimate 
\\V\n{x)-{V ~iA')ip\\\L2. 
Notice that since is radial, 

n{x) ■ V^he = 0, 

hence 

n(x) ■ (V - «A')(^ = n{x) ■ (V - iA)ip . 

Now, one deduces from Corollary 16.31 together with the energy estimate of Theo- 
rem [Ol 

(6.4) \/eP||L-(f2) < m^, ^/£\\A\\Hl(Q:) < We, \/£|| V^||l2(o) < me, 

and consequently, we obtain 

e\\A ■ V^||i2(n) < me, e|l < ■ 

Now, 

A^s = (V - iAfip + 2iA ■ Vip + \A\'^ip 

where 

{SZ-iAfip = 2— •(V-i^)^+(V-i^)-u^(V-iA)v3 

^£ 

= 2 ■ {\/ - iA')ip + {\/ - lA) ■ uliS/ - iA)ip 

[ since is radial ] , 
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and 



L°°{n) 

Therefore, we deduce from the energy estimate of Theorem lG.ll and the G-L equation 
of ip that 

e||A(p||L2(0) < nie ■ 

Now, by the standard elliptic regularity theorem, we get 

(6.5) e\\(p\\H2{n) . 
Therefore, we obtain, 

s\\(y ~iA)ip\\Hi(n) < 1- 
As we have pointed out, this is now sufHcient to deduce the first result of the lemma. 
We prove now the second statement. Let (r, 0) be polar coordinates. Let us intro- 
duce the following function in [0, 27r[, 

We claim that 

(6.6) V6ie[0,27r[, lim|/(6l)|=0. 
Actually, by the mean value theorem, we have 

\fi9)\<2£\\dM-,e)U^^[n/2,iR+i)/2]- 

On the other hand, invoking (|6.5p and the sobolev embedding theorem, we deduce 
that 

1/(0)1 «2C,^m,, 

for a constant Cg > 0. 

Ui = 0{e), this yields dHj]). Otherwise, if e{e) > e, we have 

£ £ /~2JEe\ 
-m^ < - exp <C 1, 

e s \ e J 

which again yields (|6.6p . 

Now, since \ip\ has values in [0, 1[, ||/(^, •)ll-L2([o,2ir[) is bounded. Therefore, we 
deduce by the Lebesgue dominate convergence theorem that 

1™ ll/(^>-)l|L2([0,27r[) = 0, 

E— ►0 

which is nothing but the second statement of the lemma. □ 



6.2. Proof of Theorems [TTHIIll 

Let us notice that 

"(V-iv4)'0 VWe 



il) u, 
Thus, in the regime of Theorem 16.11 we have 

'{V-iA)i}] Vu, 



\n{x) ■ (V - iA')Lp\ 



n{x) 



<2£|jn(a:).(V-a')¥>|li2(ajv)«l 



where the last conclusion is due to Lemma [ 
Invoking Theorem 13. 2^ we deduce the formulas l|1.6p and Ijl.lip . 
Now, it remains to prove that 

(6.7) II Mw|Il2(sj^) « 1- 



But, by Lemma [631 it is easy to deduce l|6.7p since tp — u^ip, and by Lemma 
Ue converges uniformly on ON to an explicit constant .4 > 0. 
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Therefore, we get now that (|1.7p and l|1.12p hold. This finally achieves the proof of 
Theorems O&O □ 



6.3. Proof of Theorem [TT6l 

Let j = {itp, (V — iA)ip) and = {iip, (V — iA)ip), where ip = 
Integrating by parts, we get 



curl = / t{x) ■ j^- I t{x) ■ , 

N J\x\ = R+i J\x\=B.-l 

where t(x) = is the tangential vector of any circle in R^. 

Notice that = z^j- Then since the function is radial and by Theorem 13.21 

ll"e-C'llL-(9Ar) ^0 (ase^O) 
for an explicit constant C > 0, we deduce that 



t{x) ■ j 



\x\=R+e 



t{x) ■ j 


= C 


/ curljy 


(l + o(l)) ase^O 


R~£ 




J N 





But, since we have no vortices, we get from Point (4) of Proposition 

/ curl = - / curlA + o(l) (e ^ 0) . 
J N Jn 

Now by the energy estimate of Theorem 16. H we deduce 



AT 



curlA < |7V|i/2||curlA|U2(n) 
< ei/2|ln£|. 
This achieves the proof of Theorem 11.61 



□ 
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Appendix A. The canonical equation 
Let P — f3{a, d) > e~^'^ and A = A{a, d) > be the solutions of the following 



equations: 



(/3eV2'i + 1)2 
. /Je^'^ + 1 ~ 



A. 



A. 
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Putting b = e^"* + e"^'' and c ^ ^/a e^'^ - ^^V^d ^ ^-^e solution (/?, A) can be 
expressed explicitly by, 

(A.l) P{a,d) = 

(A. 2) A(a,d) = 

With this choice, it is easily checked that the function U given by (|2.2p is a solution 
of Eq. (|2.ip . The aim of this appendix is to show that l|2.2p is the only bounded 
and positive solution of l|2.ip . 




Theorem A.l. In the class of functions C = {u G Hl^c{M?) n L°°{S?) : u > 0}, 
Eq. ^2.1\) admits a unique non-trivial solution 



]R2 9 {xi,X2) ^ U{X2) 

given by ^2.2]) . 



Proof. Since the proof is very close to that of [19], Theorem 1.5], we sketch only 
the main steps. 

By adjusting the proof of [191 Lemma 4.2], we obtain that if u ^ solves l|2.ip . 
then < u < 1 in R^. This permits us, when following step by step the proof of 
[19l Lemma 4.3] and [SJ, Lemma 5.3], to get a positive constant C e]0, 1[ such that 
for any solution u of l|2.ip in C, we have 



(A.3) inf uix) > C. 

Also, we prove in [191 Lemma 4.4] that, for u G C a solution of (|2.ip . 



(A.4) lim sup(l-u(a;i,a;2)) =0. 

Now, let wi, M2 € C be solutions of (|2.ip . We introduce 

(A. 5) A* = sup{A e [0, 1[ : U2(x) > Aui(x)} . 

Then, by l|A.3p . A* > 0. We claim that A* = 1. Once this is shown to hold, 

Theorem lA.il is proved. 

We argue by contradiction: If A* < 1, then 

(A.6) inf wix) =0, 

where w(x) — U2{x) — A*mi(x). Now, let (x„) — ((a;^,a;^)) be a minimizing se- 
quence: 

lim wi^Xn) = . 

n — >+oo 

Since the maximum principle yields that w{x) > for all x, the sequence (a:„) 
should be unbounded, hence we assume that lim„^+oo \xn\ = +oo. Also, by (|A.4p . 
(x^) should be bounded, hence we assume that limn^+oo x'^ = b. 
Now, the functions u"{xi,X2) = Uj{xi -|-a;", X2), j = 1, 2, solve l|2.ip in C, and up to 
extraction of a subsequence, they converge locally to functions 
€ C'ioc(K X {M \ [-d,d]}), j = 1,2. Now, ui, U2 solve dHj) in C, U2 > \*ui 
and U2(0,5) — A*ui(0,6). On the other hand, the strong maximum principle in- 
sures that U2{x) > X■^,Ul{x) for all x S M^, hence we have a contradiction. □ 
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Appendix B. On the energy of a vortex configuration with pinned 

vortices 

In this section we construct a configuration with several vortices on the circle S]^, 
and we discuss the difficulty behind the estimation of its energy. The motivation 
to construct such a configuration is that we expect vortices of a minimizer of (|1.3p 
to be pinned on the circle S]^ — dD(0,R), and to be uniformly distributed along 
this circle. 

We cover the unit disc Q — D{0, 1) uniformly by sectors, each having a small 
opening. Then we define a periodic induced magnetic field h on the sectors. This 
magnetic field will be the natural one corresponding to a vortex configuration with 
a single vortex in each sector. We mention also that similar constructions have 
been also introduced in the papers [2^i3^3|- 

Let (r, 9) be the polar coordinates. Let us take n = n{e) equidistant points 
{o-i)iLi on the circle S]^ with 9{ai) = 0. Upon defining the angular sector 

Ci = l{r,e) : re [0,1[, \e\ < 



n(e), 

we get a decomposition of £'(0, 1) by a family of disjoint sectors {Ci), where each 
Ci corresponds to the point and is congruent to Ci . 

The test configuration. 

Now we define a measure fi by: 

f if x<^U^B{a„e) 

n(x) = < 2 

^ ^ I ^ if a; e U^B{a^,e), 



and a function h' in Q — D{0,1) by 
(B.l) 

We notice that 



-div ( ^V/i' ] +h' = ^ in fl, 
h' ^0 on dn. 



^da; = 27r, V i = l,2,---N. 

c. 

We define an induced magnetic field h = h' + (here has been introduced in 
(|4.4p ). Then we define an induced magnetic potential A = A' + :^\/-^h^ by taking 
simply 

curM' = h'. 

This choice is always possible as one can take A' — V^g with g € H'^{n) such that 
Ag = h'. 

We turn now to define an order parameter if) which we take in the form 



if x G UiB(ai,e), 
if X ^ UiB{ai,2e), 

-1 if 3i s.t. X e B(ai,2e)\ B(ai,e). 

e 

The phase (j) is defined (modulo 27r) by the relation: 

(B.4) A' = -^\/^h' mn\(JiB{a„e), 

and we emphasize here that we do not need to define (j) in regions where p vanishes. 
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The energy of the test configuration. 

To estimate the energy of the test configuration {tp,A) constructed above, we ex- 
press h' by means of a Green's function, i.e. a fundamental solution of the dif- 
ferential operator —div ( V ) 4- 1. The existence and the properties of this 



(B.5) 



—div 



function, taken from |T| 128). are given in the next lemma. 

Lemma B.l. For every y G and e e]0, 1], there exists a symmetric function 
fl X n B {x,y) G^{x, y) G M+ such that : 

(^V,Ge(a;, y)^ + G,{x, y) = 5y{x) in 2?'(f7), 

Moreover, satisfies the following properties: 

(1) There exists a constant > such that 

Q<G,{x,y)<C,{\\n\x~y\\ + l), V (x, y) e f7 x 17 \ A, 

where A denotes the diagonal in IR^. 

(2) For any compact set K C il, there exist constants C > and £o > such 
that, V £ g]0, £o], 



G,{x,y) + ^^\n\x~y\ < C ' 



, y yeK.y xen. 

L°°(0) 



u1{x) 

The field h' can be expressed by means of the function G^, 

(B.6) h'{x)^ ( G,{x,y)ti{y)Ay, y x e fl. 

Jn 

By this expression of the field, one is able in the former literature to account for 
the energy resulting from the interaction of different vortices. 
By applying Lemma [6^ one essentially needs to estimate 

(B.7) f f^\s/h'f + \h'\A dx= f G,{x,y)^i{x)^l{y)dxdy. 

Unfortunately, due to the rapid oscillation of u^, the decay of G^ in Lemma [RD 
is not useful to simulate the energy, since the term || ^'^"^^ is typically of order 



Appendix C. Remarks on the Ginzburg-Landau functional 



We recall in this appendix some facts taken from [251126) concerning the Ginzburg- 
Landau functional : 

(C.l) J{u,A)^J (^\{\/ -I A)u\^ + ^{1- \u\^f + \cur\ A -H]"^^ dx, 

defined for functions {u,A) e H^{fl;C) x iJ^(f2;M^). Here 17 is a two-dimensional 
simply connected domain. 



Let us start by the following result concerning an estimate of the radius of the set 
where the modulus of a function u is small, provided that u satisfies some energy 
bound. 

Proposition C.l. There exist constants ao{fl) > and G > such that, for 
any M > 0, S > 0, e e]0, 1[ satisfying eM/d'^ < ao, any u e G^{fl;C) and 
A G G°(f);R2) satisfying 
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one has the estimate 

r{{x&Ti : \u{x)\ < 1-^}) <C^. 

The next proposition provides a lower bound of the functional J{u, A) by means 

of the degree of u on perforated sets. 

Proposition C.2. Let w be a compact subset offl. Then for any 1 > ij > a > r(w), 
there exists a collection of disjoint open balls {B{ai,ri))i such that 

(1) Ein<??; 

(2) wC[jiB{a,,ri); 

(3) For any -functions u : Cl\w and A : Q M."^ , we have 

[ \{V - iA)u\^ dx + n [ IcurlA-ifpda; > 27r|di| (in^ - ^) , 

JB{ai,ri)\w JB{ai,ri) \ a 2 J 

where di = deg{u, dB{ai,ri)) if Bi C fl and di = otherwise. 

Finally we state a result concerning 'Jacobians'. 

Proposition C.3. Assume that M > 0, B > 0, 7 e]0, 1] and S e]0,l/2]. Let 
u G C^(0;C), A e (7^(0; R^) and {B{ai,ri))i a collection of balls such that 

(1) \u\ <l inVt; 

(2) J{u,A)<M; 

(3) {x&n : < l-(5} cUiS(a^,n); 
(4) ^ri<i?. 

i 

Then, there exist measures a € H~^{rt) and (3 G (Cq (^^))' such that 
\\a\\^-^<CM5\ W\\^c^^y<CMR\ 

and, upon putting di = deg{u,dB{ai,ri)) if B{ai,ri) C and di = otherwise, we 
have, 

An ^ di - curl [{iu, V au) + A] = a + (3. 

i 

Here C > is a constant independent of e, 6, R and M. 
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